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Abstract. A method is presented for calculating the frequencies of non-retarded surface plasmons propagat- 
ing on a semi-inifinite medium with a surface profile described by a one-dimension quasiperiodic function. 
The profiles are generated, in analogy with previous work on quasiperiodic superlattices, by repeating 
unitary cells constructed according to an inflation rule. Dispersion relations are obtained for a semi-infinite 
free-electron metal as the active medium, with surface profiles obeying the Fibonacci and Thue-Morse 
sequences. 

PACS. 73.20.Mf Collective excitations - 71.45.Gm Plasmons 



1 Introduction 

Extensive experimental and theoretical work on quasiperi- 
odic systems has revealed interesting new properties which 
are not found in either periodic or in disordered systems. 
Examples of quasiperiodic systems are artificial structures 
such as surfaces with quasiperiodic tilling patterns pQ, 
as well as multilayers constructed by juxtaposing differ- 
ent building blocks following quasiperiodic inflation rules 
|2] • These structures are generally considered as interme- 
diate states between ordered and disordered systems. It 
has been found that the energy spectra of elementary ex- 
citations (e.g. polaritons, phonons, spin waves) of such 
structures are highly fragmented, with a self-similar pat- 
tern (for a comprehensive review, see Ref.0 and refer- 
ences therein). In fact, this property is now accepted as a 
signature of a quasiperiodic system. 

Another topic that has attracted considerable atten- 
tion is the study of the propagation of surface modes on 
disordered media. Among these excitations, surface plas- 
mons, which are electronic collective modes bound to a 
dielectric-metal interface, have been some of the most in- 
tensively studied. Due to their surface nature, these modes 
are strongly influenced by the shape of the interface along 
which they propagate. In fact, the excitation of surface 
plasmon modes on randomly rough surfaces can lead to 
localization effects , and is also associated with surface 
phenomena such as the enhanced second harmonic gener- 
ation 0] . Recent studies jHUHl of surface plasmons on pe- 
riodically corrugated surfaces have indicated the presence 
of absolute band gaps in the frequency spectrum, which in 
turn show a dependence on the geometry of the surface. 

The aim of this work is to present a method for cal- 
culating the spectrum of excitations of a metal-dielectric 
interface with a quasiperiodic grid of parallel ridges. In 
this case, since the quasiperiodic elements are restricted 



to the surface of the medium, they may exert a significant 
influence on its spectrum of surface excitations. Specifi- 
cally, we calculate the surface plasmon spectrum of sur- 
faces with one-dimensional profiles described by functions 
created by substitutional sequences. We use these func- 
tions to construct a series of periodic profiles, with each 
consisting of a unitary cell obtained from a quasiperiodic 
sequence, repeated along the x\ direction. For each gener- 
ation of the sequence, the period length a increases, and 
the quasiperiodic case is obtained in the limit a — > oo. 



2 Model 

The theory is based on the Rayleigh hypothesis formalism 
described in Ref. 7 , which has been shown to give ex- 
act results for surfaces with small corrugations (i.e. with 
height to width ratios < 0.072) defined by analytic func- 
tions. In this formalism, solutions of Laplace's equation 
are found that vanish as the distance (|x3|) to the inter- 
face increases. For a semi-infinite medium with a surface 
described by a one-dimensional periodic profile, the po- 
tentials associated with modes propagating along the X\ 
direction can be assumed to obey the Bloch property 



(xi + a, (J) = <p(xi,w)e 



ika 



(1) 



where a is the period of the surface. Thus, the potential 
associated with a surface mode can be expressed as, 



(xi,x 3 ) = ^C n exp(ifc„xi +a, 



X3) 



(2) 



where a n — \k + 2irn/a\, k n — k + 2irn/a, k is the wavevec- 
tor, X3 is the direction normal to the surface of the active 
medium and the C n factors are the Fourier amplitudes of 
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the potential. By applying the interface boundary condi- 
tions, the frequencies of the non-retarded surface plasmons 
can be obtained by solving the eigenvalue equation 

M rp (k)C p = — r-v — ^C r , (3) 

p—~oo 

with the elements of the matrix M rp (k) being given by 

, r /,\ a r (k)a v (k) — k r k v ,,. , , . 

M rp (k) = rv > pK ' f^-K(a r -a p ,k r - k p ) 

a r {a r -a p ) 

(4) 

for r ^ p and M rp (k) — otherwise. The function K(a r — 
a p ,k r — k p ) is obtained from the integral 

1 ,a/2 

K{a,h) = - e a ^ x) e~ lklX dx, (5) 

a J -a/2 

where we used a = a r — a p and fc; = k r — fc p . The function 
£(x) describes the surface profile. 

In the following calculations, we employ a surface model 
designed in a similar fashion to the structures used in the 
multilayer calculations. The goal is to construct a pro- 
file that approximates a quasiperiodic surface by using a 
known function, taken in finite intervals, in order to de- 
fine a periodic function. Thus a series of periodic surfaces 
can be described, each consisting of unitary cells made 
of building blocks arranged in a finite quasiperiodic se- 
quence, repeated along the Xi direction. The overall pe- 
riodic nature of the surface allows us to use the formal- 
ism of Glass et at. However, for each successive genera- 
tion of the sequence, the period length a increases, and 
an actual quasiperiodic surface is obtained as the length 
of the unitary cell grows to infinity. In the present case, 
the building blocks are chosen to be sinusoidal functions. 
These functions are defined in intervals such that each 
block corresponds to a ridge on the surface of the metal. 
This approach to the definition of the profile is similar to 
the one used in Ref [§■ for isolated ridges on a flat surface. 
The ridges can be of two types, with different heights or 
widths, which we henceforth refer to as blocks 'A' and 'B' 
(see Fig. 1). These sinusoidal blocks are described by the 
functions 

Cq (xi ) = 2A a cos 2 (ttxi I La) (6) 

in intervals s a L a + s Lb- L a /2 < x < s a L a + SbLb+L a /2, 
for 'A' ridges, and 

C b (xi) = 2A b cos 2 (7rx 1 /L b ), (7) 

in intervals s a L a + s L — L /2 < x < s a L a + SgLb + 
Lb/2, for the 'B' ridges, with s a j, = 0,1,2,.... In order 
to construct the unitary cell, one has only to specify the 
intervals corresponding to each block. 

By following the prescriptions above, one can then ob- 
tain the kernel in Eq.(3), by applying suitable changes of 
variables in Eq.(5), as 

K(a, h) = f a (h)K a (a, h) + f b {h)K b {a, h), (8) 



A 




A B A A B 




Fig. 1. Schematic representation of the surface profile used. 
Ridges, labeled A and B are arranged periodically on a flat sur- 
face. The unitary cell is constructed according to an inflation 
rule. The lower image represents a unitary cell corresponding 
to the 5th term in the Fibonacci sequence. 

where 

K c (a, fe) = - E -T / £(xi)e- ik >^dx u (9) 

with c = a, b. The functions f a {ki) and fb(ki) contain 
the information concerning the positions of the ridges in 
the unitary cell, and therefore depend on the sequence 
generation. In this paper we consider the case where A a ^ 
A b and L a = L b . 

Let us first consider the Fibonacci sequence. In this 
case, we assume So = B, S\ = A as the first two terms in 
the sequence. The other terms are generated by an oper- 
ation that consists of juxtaposing the two previous terms 
in the sequence, following the rule 

S n +i = S n S n -%. (10) 

Thus, the next three terms are 5*2 = AB, S3 = ABA, 
S4 = ABAAB etc. Therefore, for the S3 generation we 
obtain an unitary cell of length a = 2>L. This cell can be 
defined as 

rCa(xi), -3L/2 < x < -L/2, 
C(*i) = j Cb(xi), -L/2 < x < L/2, (11) 
[ Ca(xi), L/2 < x < 3L/2 

In this particular case, there is a B ridge centered at 
the origin, a A ridge shifted L units to the left of the origin, 
and another A ridge shifted by L units to the right. That 
translates as 

fb(ki) = l , /«(fcj) = 2 cos(fc,i). (12) 

For the sinusoidal profiles considered here, it can be shown 
that, 

= E|t(t) f ~" x 

n=0 



J. Milton Pereira Jr., G. A. Farias, and R. N. Costa Filho: Surface Plasmons on a Quasiperiodic Grating 3 



/2n\ sm[n(n-m)-n(r-p)/F] 
^-^Km) ir[F(n — m) — (r — p)] 

where F is the total number of ridges in the unitary cell. 

This formalism can be extended to profiles generated 
by other sequences in a straightforward way. In the case of 
the Thue-Morse sequence, one applies the following sub- 
stitution rules: 

S n = S„_r (14) 

and 

&n = ^n-l ^n— 1 (15) 

for n > 1, with So = A and Sq = B. Thus, the next three 
generations of this sequence are: S\ — AB, S2 = ABBA 
and S3 = ABBABAAB. Thus, the number of blocks in a 
unitary cell corresponding to a sequence generation S n is 
2™. 



3 Numerical Results 
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We now present numerical results for gratings ruled on a 
free-electron metal, with the dielectric constant 

e(co) = 1-%, (16) 

where lo v is the plasma frequency of the conduction elec- 
trons in the metal. The surface plasmon spectra for si- 
nusoidal gratings are known to display an infinite set of 
dispersion branches distributed symmetrically around the 
flat-surface-plasmons frequency uif sp — u> p /y2 [Jj. This 
behaviour is also found in the present results. The graphs 
in Fig. 2 show a comparison of the dispersion relations 
of surface plasmons for gratings constructed with the Fi- 
bonacci sequence, for increasing generation numbers, us- 
ing A a /L = 0.04 and Ab/L — 0.07. Only higher frequency 
branches are shown. The frequencies are normalised by 
LOf S p and the wavevectors are normalised by 2ir/a. As the 
number of blocks in the unitary cells increases, the conver- 
gence of the results becomes slower, and more terms must 
be added to the summation in Eq. (3). However, even 
for Sg, a good convergence can be obtained for N < 100 
terms. Graph (a) shows the periodic case, which corre- 
sponds to a sinusoidal grating with A/L = 0.07. Graphs 
(b) to (d) show the results for sequence generations S4 
(F — 5), S 6 (F = 13) and S 8 (F = 34), respectively. 
In contrast with the periodic case, in which most of the 
branches are found in close proximity to Wf sp , as the gen- 
eration number increases, several high-frequency branches 
are found and one can observe a tendency to the forma- 
tion of frequency bands for large generation numbers. The 
results also indicate the presence of several frequency gaps 
arising in the plasmon spectrum. This result is consistent 
with what has been previously observed for excitations 
in multilayer systems. Nevertheless, in the present case 
the width of the gaps is directly related to the height to 
width ratios of the different blocks, a parameter that can, 
in principle, be varied continuously, whereas in multilayer 



Fig. 2. Surface plasmon frequencies as a function of wavevec- 
tor for a sinusoidal grating (a), with A/L = 0.07, and gratings 
with unitary cells obtained from the Fibonacci rule: (b) S4, (c) 
S 6 and (d) S 9 , for A a /L = 0.04 and A b /L = 0.07. 

systems the gaps depend on the different properties of 
each block, which cannot be as easily controlled. 

Two results for gratings following the Fibonacci se- 
quence are displayed in Fig. 3. Both graphs show some 
of the high-frequency dispersion branches of surface plas- 
mons on gratings with unitary cells corresponding to the 
Sg term in the sequence (F = 34). However, for the graph 
labeled (a), the values A a /L = 0.04, A b /L = 0.07 were 
used, whereas for graph (b) the parameters were A a /L — 
0.07, Ab/L — 0.04. One can see from the graphs that, in 
contrast with the three high-frequency bands that arise 
in the A a > Ab case, the results for A a < Ab show the 
appearance of a single high-frequency band. These results 
illustrate the inherent asymetry of the Fibonacci sequence, 
which is a consequence of the fact that, in that sequence, 
the number of B blocks is found to be always smaller than 
the number of A blocks (for large generation numbers, the 
ratio of the number of B blocks and the number A blocks 
approaches the number w 0.618). In addition, in the 
Fibonacci sequence, the B blocks are always located be- 
tween A blocks. This asymmetry does not occur in the 
Thue-Morse sequence, in which the number of A and B 
blocks is always the same. 

Figure 4 shows the upper dispersion branches for sur- 
face plasmons propagating on a grating with a unitary cell 
corresponding to the 5th generation of the Thue-Morse se- 
quence, for A a /L = 0.02 and Ab/L = 0.07. As in the previ- 
ous cases, the graph shows several gaps in the spectrum, 
with three groups of high-frequency branches separated 
by large gaps, with several small gaps located in the high- 
frequency regions. This behavior points to a fractal nature 
of the spectrum as the generation number increases. 
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4 Conclusions 
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Fig. 3. High-frequency surface plasmon dispersion branches 
for gratings with a unitary coll corresponding to the 9th 
generation of the Fibonacci sequence, for (a) A a /L — 0.04, 
A b /L = 0.07 and (b) A a /L = 0.07, A b /L = 0.04. 
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Fig. 4. Surface plasmon dispersion branches for a grating with 
a unitary cell obtained from the Thue- Morse sequence (Sb), for 
A a /L = 0.02, Ab/L = 0.07. 



In summary, we calculated the frequencies of surface plas- 
mons propagating on gratings with one-dimensional sur- 
face profiles that follow quasiperiodic sequences. The sur- 
face morphology is described by a model that allows us 
to employ methods previously developed for the study of 
surfaces with ordered profiles, and at the same time to ex- 
plore a topography that displays deterministic disorder. 
The results were obtained for gratings constructed with 
Fibonacci and Thuc-Morse sequences and can be easily 
extended to other sequences. The surface plasmon disper- 
sions show a large number of gaps that are related to the 
sequence generation and to the aspect ratio of the ridges 
that form the grating. The spectra have a fractal character 
that is similar to what has been previously observed for 
other excitations in quasiperiodic structures. The method 
presented here allows the calculation of the frequencies 
of propagating modes. However, one cannot rule out the 
possibility of existence of non-propagating plasmon modes 
in such structures. In fact, such localized modes are often 
found in quasiperiodic supcrlattice structures, with fre- 
quencies located in the gaps of the spectrum. This method 
can also be applied to the study of other surface excita- 
tions, such as surface polaritons, and to investigate the 
light scattering properties of quasiperiodic surfaces. Sur- 
faces with periodic features have been constructed by pho- 
tolithography, and such a technique could in principle be 
applied in order to create the quasiperiodic structures de- 
scribed here. 
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